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Abstract 

We report measurements of charmonia produced in two-photon collisions and decaying to four- 
meson final states, where the meson is either a charged pion or a charged kaon. The analysis is based 
on a 395 fb _1 data sample accumulated with the Belle detector at the KEKB electron-positron 
collider. We observe signals for the three C-even charmonia g c (lS), x c o(lP) and Xc2(l-P) in the 
7r + 7r~7r + 7r~, K + K~it + it~ and K + K~K + K~ decay modes. No clear signals for g c (2S) production 
are found in these decay modes. We have also studied resonant structures in charmonium decays to 
two-body intermediate meson resonances. We report the products of the two-photon decay width 
and the branching fractions, r 77 £>, for each of the charmonium decay modes. 

PACS numbers: 13.25.Gv, 13.66.Bc, 14.40. Gx 
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I. INTRODUCTION 



The two-photon decay widths (T 77 ) of charge-conjugation (C)-even charmonium states 
are important observables that are sensitive to the properties of cc quarks inside charmonium 
bound states. Precise measurements of widths can give valuable constraints on the models 
that describe the nature of heavy quarkonia [1-3]. 

High luminosity electron-positron colliders are well suited to measurements of the two- 
photon production of charmonium states, since they provide a large flux of quasi-real photons 
colliding at two-photon center-of-mass energies covering a wide range. Various charmonium 
states, including the 7] C (1S), Xco, Xc2, Vc(2S) and x c2 (2P), have been observed so far in 
two-photon production processes. 

In resonance formation, the production cross section from quasi-real two-photon collisions 
at an e + e" collider is proportional to the product of the two-photon decay width of the 
resonance and its branching fraction (B) to the observed final state: 

a{e + e~ — > (e + e~)R, R — > final state) oc T 11 (R)B(R — > final state), 

where R is the resonance. For the above charmonium states, the two-photon decay widths 
are mainly determined from two-photon collision measurements performed in a small num- 
ber of specific decay modes only since the production rate of the corresponding states is 
limited. The results have rather large statistical errors and in some cases even larger sys- 
tematic errors originating from uncertainties in the corresponding charmonium branching 
fraction to the observed final state. It is necessary to reduce these errors with high-statistics 
measurements in various decay channels. Through these measurements, we can also obtain 
ratios of branching fractions for various charmonium decay modes. 

In this paper, we report the measurements of production of the charmonium states t) c (1S), 
Xco and Xc2 in two-photon collisions from their decays to four-meson final states 7r + 7r~7r + 7r~, 
K + K~7t + 7t~ and K + K~K + K~ , with the Belle detector. We have studied intermediate 
meson-resonant structures in these decay modes. We have also searched for production of 
the rj c (2S) in these final states. Previous measurements of the same processes are found 
in Refs. [4-10]. Belle previously reported measurements of two-photon production of these 
charmonia in several decay modes: Xc2 — > lJ/ip fl"D[, Y ^n and Xc2 decays to the two-meson 



final states 7r + 7r and K + K [12j as well as KgK$ 13], and rj c (lS) — > pp [14|. Production 



of the T] C (2S) in a two-photon process was observed in the KgK^ir mode by the CLEO 



and BaBar experiments [15j, [16j. This is the only decay mode so far directly observed for 
this charmonium state. 



II. EXPERIMENT AND EVENT SELECTION 
A. Data and the Belle detector 

We use data that corresponds to an integrated luminosity of 395 fb _1 recorded with 
the Belle detector at the KEKB asymmetric-energy e + e~ collider [17] . Since the beam 
energy dependence of two-photon processes is very small, we combine the on-resonance and 
off-resonance data samples; the off- resonance data were taken 60 MeV below the T(45) 
(-^/s = 10.58 GeV). The analysis is made in the "zero-tag" mode, where neither the recoil 
electron nor positron is detected. We restrict the virtuality of the incident photons to be 
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small by imposing a strict transverse-momentum balance along the beam axis for the final- 
state hadronic system. 



A comprehensive description of the Belle detector is given elsewhere [18j. We mention 
here only the detector components essential to the present measurement. Charged tracks are 
reconstructed from hit information in a central drift chamber (CDC) located in a uniform 
1.5 T solenoidal magnetic field. The z axis of the detector and the solenoid are along the 
positron beam direction, with the positrons moving in the —z direction. The CDC measures 
the longitudinal and transverse momentum components (along the z axis and in the rip plane, 
respectively). Track trajectory coordinates near the collision point are provided by a silicon 
vertex detector (SVD). Photon detection and energy measurements are performed with a 
CsI(Tl) electromagnetic calorimeter (ECL). Kaons are identified using information from 
the time-of-flight counters (TOF) and silica-aerogel Cherenkov counters (ACC). The ACC 
provides good separation between kaons and pions or muons at momenta above 1.2 GeV/c. 
The TOF system consists of a barrel of 128 plastic scintillation counters, and is effective in 
K/tt separation mainly for tracks with momentum below 1.2 GeV/c. Lower energy kaons 
are also identified using specific ionization (dE/dx) measurements in the CDC. The magnet 
return yoke is instrumented to form the Kl and muon detector (KLM), which detects muon 
tracks and provides trigger signals. 

Signal events are efficiently triggered by several kinds of track-triggers that require two 
or more CDC tracks with combinations of TOF hits, ECL clusters or summed energies 
of ECL clusters. No additional cuts on the trigger information are applied. The trigger 
conditions are complementary for the detection of four-prong events; we obtain an overall 
trigger efficiency of ~ 95%. 

Kj 7r separation uses a likelihood ratio formed from ACC, TOF and CDC information. No 
explicit lepton identification requirement is applied, since leptons are not a large background 
source in exclusive four-prong events. 



B. Event selection 

Candidate events are selected as follows; all variables in selection criteria (1)-(10) are 
measured in the laboratory frame. (1) There are exactly four charged particles and each 
satisfies the following criteria: pt > 0.1 GeV/c, dr < 5 cm, \dz\ < 5 cm, where pt is the 
transverse momentum of a track with respect to the z axis, and dr and dz are the radial 
and axial distances, respectively, of the closest approach (as seen in the rip plane) to the 
nominal collision point; (2) the net charge of the four tracks is zero; (3) at least two of 
the four tracks satisfy the following additional criteria: p t > 0.4 GeV/c, dr < 1 cm, and 
—0.8660 < cos# < +0.9563, where 9 is the polar angle; in addition, the following two criteria 
are applied for the four-track systems: (4) the scalar sum of the momenta \pi\) of the 
four tracks must be smaller than 6 GeV/c; (5) the average of the dz values for the four 
tracks is within 3 cm of the nominal collision point. 

In order to reject backgrounds from single-photon annihilation or radiative events, the 
following criteria are required: (6) the invariant mass of track combinations that satisfy 
criterion (3) should be smaller than 4.5 GeV/c 2 (here, zero mass is assigned to each of the 
charged tracks), and (7) the missing mass squared of the recoil against the track combination 
is larger than 2 (GeV/c 2 ) 2 . 

Calorimetry requirements are also applied: (8) the total calorimeter energy in an event 
is smaller than 6 GeV; (9) there is no electromagnetic cluster unassociated with a track 



5 



x 10 3 



3500 

_o 

> 3000 

^ 2500 

° 2000 

| 1500 

§J 1000 
LU 

500 







0.05 0.1 0.15 

|Zptl (GeV/c) 



FIG. 1: Distribution of the vector sum of transverse momenta. Events to the right of the arrow 
are discarded. 

(i.e., a neutral cluster) whose energy exceeds 0.4 GeV; (10) there is no it candidate whose 
transverse momentum is larger than 0.1 GeV/c. 

After all these selection criteria are applied, the four tracks are transformed into the 
e + e~ center-of-mass (cm.) frame, and the vector sum of their transverse momenta with 
respect to the beam direction in the cm. frame, £pj\ is calculated. The variable Ep^ 
approximates the transverse momentum of the two-photon-collision system. We require 
(11) | Sp* | < 0.1 GeV/c in order to enhance the number of events from quasi-real two- 
photon collisions. Figure 1 shows the transverse momentum (|£p£|) distribution near the 
selection region. 



C. K/tt separation 

We apply particle identification requirements for each track. A likelihood ratio, 1Z = 
Lk 1 1 (Lk+L^), is calculated, where the kaon and pion likelihoods, Lk and L n , are determined 
from information provided by the ACC, TOF and the CDC dE/dx system. All tracks with 
1Z > 0.8 are assigned to be kaons. This requirement gives a typical identification efficiency 
of 90% with a probability of 3% for a pion to be misidentified as a kaon. All other tracks 
are treated as pions. 

We take (12) only the combinations in which the net strangeness is conserved, 7r + 7r~7r + 7r~, 
K + K~7r + ir~ and K + K~ K + K~ . These final-state combinations are sometimes denoted 4tt, 
IKI-r and 4K, respectively, in this paper. 

When a pair of pions (kaons) in the four prong decay of a certain charmonium are 
misidentified as a pair of kaons (pions), they could produce a broad enhancement in the 
spectrum of the misidentified channel at a mass very different from that of the original 
charmonium. However, these backgrounds do not make any pronounced peaks near the 
three charmonia masses. 
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FIG. 2: A two-dimensional plot for T,p z vs invariant mass of the hadronic system for the four-pion 
final state candidate. The curve is the cut criterion (13) described in the text. 

D. Rejection of ISR events 

Two-photon events with relatively high 77 cm. energies (above 3.2 GeV in the case of 
Belle) are contaminated by background processes where a real or virtual photon is emitted 
in the direction of the incident electron. They are called ISR events or pseudo-Compton 
events, in case the positron collides with the electron (e + e~ — > (e + e~*)7 — > hadrons + 7) or 
with the virtual photon (e + e~ — ► (e + 7*)e~ — > hadrons + e + e~), respectively. Such events 
have a kinematical correlation between the p 2 -component in the laboratory frame (Ep z ) and 
the invariant mass of the final-state hadron system iW). We reject them with the following 
additional requirement: (13) Ep z > (W 2 - 49 GeV 2 /c 4 )/(14 GeV/c 3 ) + 0.6 GeV/c. These 
backgrounds are not harmful in the measurements of C-even charmonia, since at lowest 
order they produce only C-odd hadron systems. We show a two-dimensional plot for Ep 2 
vs invariant mass of the hadronic system for the four-pion final state candidate in Fig. 2. 

After all these selections, we obtain 6.37 x 10 6 47r events, 4.17 x 10 5 2K2ir events and 6003 
4K events in the four-meson invariant mass region between 2.0 GeV/c 2 and 5.0 GeV/c 2 . 
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III. YIELDS FOR CHARMONIUM PRODUCTION 



A. Yields of the charmonia 

We measure the charmonium yields from the invariant mass distributions in each of 
the three final-state processes. The distributions are shown in Fig. 3. Peaks at the C- 
even charmonium masses are assumed to be from two-photon production processes if the 
transverse-momentum distribution peaks in the vicinity of |£p£| = 0. We discuss the |£p£| 
distributions in Sect. V.C. 

We find clear enhancements for the rj c (lS) (denoted hereafter as "?? c ") at ~ 2.98 GeV/c 2 , 
Xco at ~ 3.41 GeV/c 2 , and Xc2 at ~ 3.555 GeV/c 2 in all the final states. We do not find 
any clear rj c (2S) signals, which would appear at around 3.63-3.65 GeV/c 2 in the invariant 
mass distributions. A small peak near 3.69 GeV/c 2 in the Air final state is from ip{2S) — > 
J/ijj7r + 7r~ , J/iJ; — > where the leptons are not identified and are treated as pions. 

These events are attributed primarily to the double ISR e + e~ annihilation process, e + e~ — ► 
i/)(2S)yy. Upper limits for rj c (2S) production will be derived in Sect. VII. 

The invariant-mass distribution in the vicinity of each charmonium peak is fitted to the 
sum of charmonium and background components using \ 2 minimization. We use a second- 
order polynomial for the background component. The fitted regions are 2.8-3.2 GeV/c 2 , 
3.3-3.5 GeV/c 2 and 3.5-3.6 GeV/c 2 for the analyses of the T] c , Xco and Xc2, respectively. The 
energy dependence of the efficiency is small within each resonance region and is neglected. 
We adopt the relativistic Breit-Wigner function for the r\ c and Xco signals, which are smeared 
with a Gaussian function whose width is fixed to the invariant-mass resolution (ranges in 
7-10 MeV/c 2 depending on the decay process) estimated by the signal Monte Carlo (MC) 
simulation. In the fit to the Xc2, we use a simple Gaussian function with a floating width, 
because the natural width is considerably smaller than the mass resolution. The mass 
resolutions from the fits, 8-9 MeV/c 2 , are consistent with the MC expectations at the Xc2 
mass. 

The fit results are shown in Fig. 4 and are summarized in Table [H In the table, the 
masses are corrected for the effects of a systematic shift seen in the signal MC mainly due 
to the tails of kaon energy loss in the innermost detector region. The correction size is 
typically x (1 - 2MeV/c 2 ), where nx is the number of kaons. The systematic errors for 
the masses include the uncertainty due to this effect (a half of the correction size) as well 
as the uncertainty of the mass scale (2.0 MeV/c 2 ) originating from the uncertainty of the 
momentum scale, while those for the widths are obtained by changing the mass resolution 
by 1 MeV/c 2 . We confirm that the mass and width of the 7] c are robust for the change of 
the fit region as described in Sect. V.B. 

We find a hint of a possible J /if) contribution from the double ISR background in the 
2K2-R distribution. This could affect the yield and resonance parameters of the r\ c determined 
from the fit. We have tried a fit that includes a peak near the J /if} mass. The changes of 
the i] c yield (7%) and the total width (1.9 MeV) are taken into account in the systematic 
errors. There is no detectable shift in the mass. The fit including a J /if) peak gives the J / if) 
mass, 3093.6 ± 2.0 ± 2.4 MeV/c 2 , where we apply the same correction and systematic error 
as applied for r) c —* 2K2n. This provides a confirmation of the mass scale in the present 
measurement. 

Possible interference of the charmonia with the continuum component is not taken into ac- 
count in the present analysis. It is difficult to separate the continuum-component amplitudes 
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FIG. 3: Invariant mass distributions for the four-meson final states, (a)47r, (b)2i^27r and {c)AK. 



that interfere with the charmonium amplitude in each process, because different partial- wave 
components of the continuum appear within the background. In general, the interference ef- 
fects could give different mass/width fit results for each decay mode. However, for the three 
considered processes the values obtained for these parameters from the no-interference fits 
to the separate decay modes are consistent with each other. The mass values for the rj c are 
systematically higher than the world average value, 2979.8 ± 1.2 MeV/c 2 Q. The r] c width 
and the parameters of Xco and Xc2 are consistent with the previous measurements. Table [I] 
also shows the combined results from the measurements of the three decay processes, where 
the central values are derived from the weighted means of the three fits and the systematic 
errors are evaluated considering the full correlations between the error sources. 



B. Yields of two-body decays 

We study the sub-structure of the charmonium signals by searching for quasi-two-body 
components, in which a charmonium meson decays into two resonances, and each resonance 
decays to two final-state mesons. We determine the yields of these two-body decays using 
the following procedure: we first make two-dimensional plots of the two invariant masses 
constructed from the available two-meson combinations of tc + it~ , K ±r R T or K + K~ , for each 
event. The An and 4K samples have two entries from each event. Such plots are made in 
each charmonium signal region and sideband region. The charmonium sideband distribution 
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FIG. 4: Results of the fits to the invariant mass distributions for (a,b,c) 4-7T, (d,e,f) 2K2ir and 
(g,h,i) AK final states in the charmonium mass regions. 



is then subtracted from the signal distribution bin-by-bin in two dimensions. We note that 
the non-charmonium background components are considerable even in the vicinities of the 
charmonium peaks. 

The signal regions are defined as the ranges within ±50 MeV/c 2 (±30 MeV/c 2 for the 
Xc2) of 2.98 GeV/c 2 , 3.41 GeV/c 2 , and 3.56 GeV/c 2 , for the r] c , Xco and Xc2, respectively. 
The sideband regions are taken on both sides of the signal region, [—100, —50] MeV/c 2 below 
and [±50, ±100] MeV/c 2 above ( [-60,-30] MeV/c 2 and [±30, ±60] MeV/c 2 for the Xc2) 
the central point of the corresponding signal region. The sizes of the signal and sideband 
regions (in the sum of the two regions) in the invariant-mass range are chosen to be the 
same so that we can apply a simple sideband subtraction. 

Next, we define the signal and sideband regions for the p° — > 7r + 7r± / 2 (1270) — ► 7r + 7r± 
i^*(892)° — ► K ± w T and <p — > K + K~ mesons, based on their known masses and widths. 
(Hereafter, we refer to the / 2 (1270), K*(m2)° and /£(1525) as f 2 , K*° and fa respectively.) 
We apply the sideband subtraction at the / 2 side to extract the yield of the / 2 / 2 decay. 

The signal regions are 0.64-0.88 GeV/c 2 , 1.08-1.40 GeV/c 2 , 0.80-0.96 GeV/c 2 and 1.00- 
1.04 GeV/c 2 for the p°, f 2 , K*° and 0, respectively. We take sideband regions on both sides 
adjacent to the signal region. In the <fr K + K~ mode the sidebands cover the ranges from 
the mass threshold to 1.00 GeV/c 2 and 1.04-1.08 GeV/c 2 . For each of the other resonances 
the sideband has a width that is half the size of the signal region on each side. The sideband 
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TABLE I: Results for the resonance parameters obtained from the fits to the invariant-mass dis- 
tributions. The third column shows the Breit-Wigner width for the r/ c and XcO> an d the Gaussian 
standard deviation for the Xc2- The combined results from the three decay processes are also shown. 
The first and second errors for the masses and widths are statistical and systematic, respectively. 
The other errors are statistical only. 



Charmonium — * 


Mass 


Width(T) (or a M ) 


Number of events 


X 2 /d.o.f 


final state 


(MeV/c 2 ) 


(MeV) 






Vc ~ 




2987.5 ±1.7 ±2.0 


36.0 ± 5.0 ± 1.7 


5381 ± 492 


48.3/34 


Vc -> 


2K2tt 


2983.6 ± 1.6 ± 2.4 


23.2 ±5.1 ±2.8 


2019 ± 248 


61.7/34 


Vc ~ 




2988.0 ±2.3 ±4.0 


19.0 ±7.8 ±2.2 


216 ± 42 


35.5/34 


XcO 


— > 47T 


3414.7 ±0.7 ±2.0 


11.1 ±2.6 ±2.3 


3550 ± 273 


27.1/14 


XcO ~ 


> 2K2ir 


3413.0 ±0.8 ±2.3 


10.1 ±3.0 ±3.1 


1694 ± 161 


24.8/14 


XcO ' 




3415.6 ±2.2 ±3.1 


10.1 ±7.1 ±2.2 


215 ± 36 


24.7/14 


Xc2 


— ► 47T 


3554.4 ±0.9 ±2.0 


(9.2 ±0.8) 


1597 ± 138 


4.7/4 


Xc2 ~ 


> 2K2-K 


3555.0 ±1.0 ±2.1 


(8.4 ±0.8) 


780 ± 74 


2.8/4 


Xc2 


-» 4K 


3558.0 ±1.8 ±2.9 


(8.1 ± 1.3) 


126 ± 24 


3.9/4 


T] c (combined) 


2986.1 ± 1.0 ±2.5 


28.1 ±3.2 ±2.2 






XcO (combined) 


3414.2 ±0.5 ±2.3 


10.6 ±1.9 ±2.6 






Xc2 (combined) 


3555.3 ±0.6 ±2.2 









subtraction is applied assuming the continuum yield has a linear dependence on invariant 
mass. 

Thus, we obtain a one-dimensional histogram corresponding to a resonance signal com- 
ponent tagged by another resonance. We search for p°, /a, K*°, <fi and f' 2 — * K + K~ signals 
in the histograms. In the decays to pairs of same-kind resonances (such as ^2/2 or K*°K*°), 
we symmetrized the two-dimensional distribution in two directions, by adding it to the 
same distribution where the two-masses are interchanged. In this case we must divide the 
peak yield by two to obtain the correct signal yield and take into account the statistical 
correlation. 

We fit the signal to a Breit-Wigner taking into account the angular momentum between 
the final state mesons and a linear continuum component in the one-dimensional distribu- 
tion near the resonance masses. As the only exception, however, we count in <p<fi decays the 
number of events in the K + K~ invariant mass range, 1.00-1.04 GeV/c 2 and subtract back- 
grounds from a sideband in the 1.04-1.24 GeV/c 2 region, instead of performing a fit. The 
one-dimensional distributions and the fits for the p°, /2, K*° and f 2 resonances in different 
charmonium decay modes are shown in Figs. 5-7. 

Table [IT] summarizes the signal yields. We assume that there are no effects from inter- 
ference in the sideband subtractions as well as the fits of the invariant mass distributions. 
We have applied corrections for inefficiencies arising from sideband subtractions in the two- 
meson resonance and charmonium regions; the targeted resonance component partially leaks 
outside the signal region and enters the sideband regions. The inefficiencies are calculated 
by assuming Breit-Wigner forms for the charmonia and light-quark resonances using the 
known masses and widths 19| . The 90%-confidence-level (C.L.) upper limits correspond 
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to the upper edge of a 90% confidence region derived from Feldman and Cousins's table 
(Table X in Ref. 20] ) assuming statistical fluctuations follow Gaussian distributions (there 
is a large number of events in both signal and sideband regions in the present case), with 
the expected number of signal events constrained to be non-negative. 

We do not observe a significant p°p° component in r] c decays in spite of the signal reported 
for r] c — > pp from measurements of radiative J/if> decays [19J. Meanwhile, we find and 
72/2 signals with greater than Aa statistical significances. It can be noted that the observed 
yield for r\ c — ► ^2/2 — ► 47r is comparable to the total yield of r\ c — > An. This indicates 
that the fifi component dominates r\ c — > An decays. We confirm the decays to jifi in 
the distribution of the transversity angle (Fig. 8(a)) defined as the angle between the decay 
planes of the two f'2 candidates in the r] c rest frame. Here, the ji candidate is a n + n~ 
combination whose invariant mass lies in the range 1.15-1.39 GeV/c 2 . The distribution 
shows a characteristic feature of tensor-meson decays. Similarly, rj c —>■ 2K2tt is dominated 
by the sum of r) c — ► K*°K*° — > 2K2n and t] c — > f 2 f% ~~ > 2K2n. The transversity-angle 
distribution of the ^2/2 candidates is shown in Fig. 8(b) (here, the f 2 candidate is a K + K~ 
combination whose invariant mass falls in 1.475 - 1.575 GeV/c 2 ). The asymmetry with 
respect to 90° can be interpreted as the effect of interference between the above two decay 
modes; an asymmetry in directional correlations between K + n + and K + n~ should not arise 
in a pure f 2 fy process. The (M(K~n + ), M(K + n~)) distributions expected from the fif'i MC 
events overlap with the (K*°, K*°) mass region near its lowest mass-combination edge. The 
histograms are the distributions of the corresponding signal MC events generated assuming 
pure P-waves between the two tensor mesons. 

Figures 9 (a) and (b) show the invariant-mass distributions of An events that are consistent 
with p°p° and ^2/2, respectively. We find no prominent peaks at the r\ c mass in the former, 
but find a clear signal for the 7] c in the latter. An enhancement at the 7] c mass is also clearly 
visible in Fig. 9(c) which shows the invariant-mass distribution for fif^ — ► 2K2n candidates. 



C. Yields of three-body decays 

The distributions corresponding to the processes Xc2 —> p°vr + 7r~ in Fig. 5(c) and 
Xco — > K*°K~tt + or K*°K + tt~ in Fig. 6(b) have significant net yields although there are no 
prominent structures from the two-body decay components, Xc2 — * P°P° and Xco — > K*°K*°. 
We obtain the yields of each three-body decay using the following procedure. For the two 
processes, we first obtain the numbers of p°'s in Xc2 —> 47r and K*°'s in Xco 2K2tt by 
fitting the invariant-mass distributions of 7r + 7r~ (Fig. 10(a), 4 entries/event) and K ± 7r T 
(Fig. 10(b), 2 entries/event) to functional forms corresponding to the sum of p° and fi res- 
onances and a K*° resonance, respectively. Here, we also take into account the continuum 
component parameterized by a second-order polynomial. The results obtained are converted 
into the yields of three-body decay events after subtracting possible resonance yields from 
Xc2 p°7i + 7T~ via p°p° and Xco — * K*°K + 7t~ (or c.c.) via K*°K*° using the best estimates 
of the corresponding two-body yields after doubling the number of events to take account 
of the multiple entries. The results after these subtractions are summarized in Table [Til 
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TABLE II: The yields of two-body and three-body decay events. The errors are statistical. The 
efficiency due to the sideband subtraction procedure is corrected. The upper limits are at the 
90% C.L. 



Decay process 


Number of events 


n n 
rj c -» p u p v - 


-> 47T 


< 1556 


Vc -> HH - 


-> 47T 


OloZ ± / DO 


n n 

XcO -» /5 U P U 


— > 47T 


^ oro 
< ZOZ 


Xc2 — * P°P° 


1 /I rrr 

— t Q7T 


^ rno 
< 0»O 




-> 2K27T 


882 ± 115 


r? c -> /2/2 -> 


2K27T 


1128 ± 206 




-> 2K27T 


< 148 


Xc2 - *° 


-> 2K27T 


151 ± 30 


77c -> 00 — 




132 ± 23 


XcO 00 - 


+ 4K 


23.6 ±9.6 


Xc2 ~> 00 - 


* 4K 


26.5 ±8.1 


Xc2 -> p°vr + vr" 


~ — > 4ir 


986 ± 578 


XcO - J ST*°^-vr+(or 


C.C.) -> 2K27T 


495 ± 182 



IV. DERIVATION OF PRODUCTS OF TWO-PHOTON DECAY WIDTH AND 
BRANCHING FRACTION 



The product of the two-photon decay width and branching fraction is obtained from the 
following formula: 

N 

Q{R — > final state) = T 7 j{R)B(R — ► final state) = — — - — — — -, 

e J Cat F(M R , J) 

where N is the number of observed events, e the efficiency, and / Cdt the integrated lu- 
minosity. F(Mr, J) is a factor that is calculated from the two-photon luminosity function 
(L 77 (Mfl)), the charmonium mass (Mr) and the charmonium spin (J), using the relation 

F(M R , J) = 4tt 2 (2J + l)L^{M R )/M R \ 

that is valid when the resonance width is small compared to the available kinetic energy 
in the decay. The F parameter is calculated to be 2.11 fb/eV, 1.15 fb/eV, 4.74 fb/eV and 
0.861 fb/eV for the T] c , Xco, Xc2 and r] c (2S), respectively, using the TREPS code [2l|. We take 
the signal yields from Tables HI and ILT1 We subtract the small contributions from Kg Kg decays 
included in the Xco ~^ 47r and Xc2 - ► 47r samples, 142±6 and 35±5 events respectively, which 
were determined from the two-dimensional invariant-mass plots mentioned in Sect. III. A. 
(XcO; Xc2 - > K S K S results from our experiment are reported separately [l3j. r\ c — > K S K S is 
prohibited from parity conservation.) 

The efficiency in the detection of each process is determined from the signal MC events 
generated by the TREPS code Eil and processed by a full detector simulation. Three- 



and four-body decays were generated according to phase space with a resonance mass dis- 
tributed according to the Breit-Wigner function, and the K*°(p°) decaying to K ±r K T (7r + 7r~) 
isotropically. 
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FIG. 5: The two- meson invariant-mass distri- 
butions and fits used to determine two-body 
yields in (a) 7] c -> p°p°, (b) XcO -> P°/A (c) 
Xc2 -> P°/0°, and (d) r? c -> / 2 / 2 . A p° or a / 2 
is selected on the opposite side. We rescale the 
number of entries described in Sect. III.B by a 
factor of 1/2 to normalize the vertical axis to 
the number of events. 



FIG. 6: The two-meson invariant-mass distri- 
butions and fits used to determine two-body 
yields in (a) r/ c -> K*°K*°, (b) X c0 -> K*°K*°, 
(c) X c2 - and (d) r, c - / 2 / 2 . A 

ET*° or a / 2 is selected on the opposite side. 
We rescale the number of entries described in 
Sect. III.B by a factor of 1/2 to normalize the 
vertical axis to the number of events, except 
for (d). 
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FIG. 7: The two-meson invariant-mass distributions in (a) rj c — ► ^0, (b) XcO — * 00) an d (c) 
Xc2 - > 00- A is selected on the opposite side. We rescale the number of entries described in 
Sect. III.B by a factor of 1/2 to normalize the vertical axis to the number of events. 



For two-body decays involving resonances with spin, the decay amplitude takes into 
account correlations caused by spin, and the full matrix element is symmetrized with respect 
to identical final-state particles. We assume a pure state with the lowest possible orbital 
angular momentum that satisfies conservation laws in each decay channel. The helicity 
of the Xc2 along the incident axis is assumed to be 2, following the previous experiment 
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transversity angle (rad) 



FIG. 8: The transversity-angle distributions for (a) rj c — > / 2 / 2 — > 47r and (b) ?? c — > /2/2 — > 2K2-K 
candidates for the experimental data after the subtraction of charmonium sideband contributions 
(closed circles with error bars). The histograms are the distributions from signal MC (see the text 
in Sect.III.B). 




M(4tt) (GeV/c2) m(4ti) (GeV/c 2 ) M(2K2ti) (GeV/c 2 ) 

FIG. 9: The An invariant-mass distributions for (a) p°p° and (b) ^2/2 candidates and the 2K2ir 
invariant-mass distributions for (c) 72/2 candidates. No subtraction is applied for p°, / 2 and 
sidebands, that is, contributions from p tt + -k~ etc. may also be included. 




M(nn) (GeV/c2) M(Kn) (GeV/c2) 

FIG. 10: (a) it + it~ invariant distributions in the Xc2 decay products and (b) K^tt^ invariant 
distributions in the XcO decay products. The charmonium-sideband distributions are subtracted. 
The curves are the fits considering the resonances decaying into two mesons (see text in Sect. III.C). 
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results 11] and theoretical expectations 22l . 

The trigger efficiency is taken into account; it is estimated to be around 95% within 
the acceptance for all the processes. The detection efficiency including that of particle 
identification is calculated using MC events passed through the detector simulation code. 
Typical efficiencies are 16%, 10%, 5%, and 9% in the 47r final-state processes, IKI-a final- 
state processes, 4K with a phase-space distribution, and the 00 process, respectively. Their 
dependence on the two-photon cm. energy, or charmonium mass, is rather small. 

We summarize the results on Q = r 77 £> in Table IHIi The systematic errors and com- 
parison to the previous experimental results also given there are described in detail in the 
following two sections. In some decay channels, we assume isospin invariance and use nec- 
essary branching fractions from Ref. [l9[. The "isospin -BF" factors in the table are the 
products of the isospin factor and the branching fraction(s) used to obtain the results. 



V. SYSTEMATIC ERRORS 
A. Uncertainties in efficiencies 

The trigger efficiency is estimated using the Monte-Carlo trigger simulator to be (95 ± 
4)%, where the error is estimated from comparison of experimental and MC results in low- 
multiplicity events and the efficiency variation under different beam background conditions. 

We obtain a 5.5% systematic error from the uncertainty in track reconstruction efficiency 
for the four-track events, as well as 2% per kaon track from the kaon identification efficiency. 
The latter is determined from a study of kinematically identified kaons in a sample of 
D* + — > D°7i + , D° — > K~-k + decays. We check the kaon identification efficiency by loosening 
the 4A"-event selection requirement; we temporarily require only three kaons in an event. We 
obtain (65±16)% more yield for the rj c , Xco and Xc2, on average with the looser condition. The 
charmonium peaks are less prominent with these selections, due to a larger contamination 
from non-4A" events. The above ratio is consistent with expectations from the MC, ~ 53%, 
for the average of the three charmonia decaying to the 4K state. 

We use the total widths of the charmonia for the determination of the efficiency in the 
sideband subtractions. Taking systematic errors in the r\ c and Xco widths into account, we 
varied the widths of the r] c and Xco from their nominal values by ±10 MeV and ±7 MeV, 
respectively, and assigned the obtained efficiency changes, ±12% and ±8%, respectively, as 
a systematic error. We take a 3% systematic error in the Xc2 case for the uncertainty in the 
treatment of its long Breit-Wigner tails. 

Inefficiencies caused by the sideband subtractions of the isobar resonances decaying to 
two mesons are calculated using the resonance parameters. We estimated the systematic 
error from this effect by changing the integration range of the resonance function from ±8r 
to ±5r. It is estimated to be 5% for each two-body decay mode. 



B. Systematic errors in signal yields 

In the fits to the four-meson invariant-mass spectra described in Sect. III. A, we set the 
fit region for the analyses of the 7] c to be rather wide, considering the relatively large natural 
width and the relatively large size of the continuum components. As a result, we have a 
larger second-order polynomial term contribution from the continuum components, which 
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correlates significantly with the size of the charmonium component and could bias the signal 
yield. We tentatively narrowed the fit region to 2.85 - 3.15 GeV/c 2 and took the variations 
in the obtained yields as a systematic error. They amount to 7%, 13% and 24% for 4ir, 
2K2-K and 4K channels, respectively. We neglect this error source for the other charmonia, 
where the quadratic components in the continuum are small. 

In the fits to the two-meson invariant-mass distributions used to obtain the yields of the 
two-body decays, we tried a second-order polynomial constrained to vanish at the two-body 
threshold for each continuum component. We assign the difference between the yields from 
this method and the standard fit (to a linear function) to the systematic error. Thus, the 
systematic errors are 20% for rj c — > f 2 f^ and less than 3% for the other two-body channels 
where we have significant signals. 

In the case of the three-body modes, we tested an alternative functional form for the 
continuum component and used the shape of three-body phase space calculated by the MC 
generator and the detector simulator. We take the difference in the signal yields from this 
method as the systematic error, that is 4% and 7% for Xc2 — > p n + n~ and Xco K^K^ir^, 
respectively. 

C. Study of non-exclusive backgrounds 

Backgrounds can come from general multi-hadron production in two-photon collisions 
and qq production in single-photon annihilation. Yields from the cc process in the latter 
category are estimated from the MC and found to be negligibly small (~ 1%) compared to 
the non-charmonium backgrounds from two-photon processes. However, the charmonium 
production rates from these processes cannot be reliably estimated. 

A more reliable signature is the pt-balance nature in the signal events. We expect that 
non-27 background processes cannot produce a p t -balanced distribution, except for double- 
ISR processes where exclusive production of C = + charmonia is prohibited. Even in the 
decays i]){2S) — > 7Xcj where a tp(2S) is produced in an ISR or double ISR process, the p t 
distribution of Xcj will not peak inside the selected region, |Ep£| < 0.1 GeV/c. We estimate 
the yield of this process to be about 1% of the Xco and Xc2 signals, based on the observed 
yields of ip(2S) decaying to 7T + 7T~J/ifj. 

We evaluate the non-exclusive backgrounds by examining the observed |£p£| distribu- 
tions. Figure 11 shows the |£p£| distributions of the charmonium signal component after 
the sideband subtractions. The distributions are compared with the signal-MC expectations 
normalized to the sum of the two leftmost bins where the background contamination is ex- 
pected to be very small. We expect that the background component has a linear shape in 
|Sp^|, in the region below 0.1 GeV/c, based on geometrical considerations and some MC 
studies. Given this assumption, we obtain the background level from data, to be (6 ± 7)%, 
(—4 ±5)% and (9 ±7)%, on average for the r) c , Xco and Xc2 decay modes, respectively, using 
the experimental excesses seen in the bins of 0.07 < |Sp*| < 0.10 GeV/c. Since there is no 
significant difference among these estimates, and they are consistent with zero (the overall 
average is (2 ±4)%), we do not apply any corrections to this background source and assign 
a systematic error of 6% for all processes. 
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FIG. 11: The | I distributions for the charmonium signal components after the charmonium- 
sideband subtractions (closed circles with error bars), compared with the signal MC (histograms) 
in the processes (a)rj c ->■ 4tt, (b)xc0 -> 47r, (c)xc2 ->• 47r, (d)r/ c -> 2K27T, (e)xco ->• 2ET27T, (f)%c2 -> 
2K2ir, (g)r/ c — > 4iC (h)xco — ► 4X and (i)xc2 - ► 4ET. The MC is normalized to the sum of the two 
leftmost bins. 



D. Total systematic errors 

We consider other sources of systematic errors. The statistics of the signal MC events 
gives errors of 2-4% depending on the process. The calculation of the luminosity function 
has a systematic uncertainty of 5%. We take the differences among the efficiency estimates 
for different spin/angular- momentum assumptions for the Xc2 decays to K*°K*° and <p(f) as 
systematic errors. The nominal S-wave hypothesis gives the most conservative upper limit 
for the Xc2 P°P° case. The total systematic errors are listed in Table III. The size of the 
systematic error relative to the central value is process dependent and ranges from 11% to 
35% in the processes where finite Q results are presented. The upper limits in the table 
are shifted upwards by the la systematic errors, in order to include the effect of systematic 
uncertainties. 



VI. COMPARISON WITH PREVIOUS MEASUREMENTS 

Some previous measurements give Q results for the processes measured here a a a. 
We make direct comparisons with these measurements, which are listed in Table IIII1 We 
quote the world average value from Ref. 19| when two or more previous measurements are 
available. We can make such direct comparison only for a limited number of processes, since 
our results include many processes measured for the first time. 

It is possible to make some indirect comparisons by converting the Q results measured 
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for a different decay mode from the present measured mode and multiplying by the ratio of 
the branching fractions; 

where R — > A is a "normalization" process for which previous measurement (s) are available. 
We adopt r\ c — > KKn, Xco -k + h~ / K + K~ (the average of the two modes) and Xc2 iJ/^P 
as the normalization processes. Note that the two-photon decay widths for these charmonium 
states are so far measured only in a few decay modes; there are not enough consistency checks 
among the different decay modes available so far. We use the PDG's fit when the average 
is unavailable and treat the errors shown there as independent from each other. 



We find large differences between the present and previous measurements for i] c in the 
direct comparisons, although they are not inconsistent considering the large errors assigned 
to the previous measurements. However, in the indirect measurements, for many r\ c decay 
modes our Q values are consistently smaller by factors of two to four compared to the previous 
determinations. We do not observe any significant signal in the p°p° mode in contrast with 



previous results |23|, ]24j. We conclude that these discrepancies are not likely to be due to 
statistical fluctuations, and there could be systematic deviations in previous measurements 
of the relevant partial decay widths or branching fractions. 

In contrast, all the present Xco and Xc2 results are in agreement with previous values, and 
confirm the values of r 77 (x c0 ) and r 77 (x c2 ) results derived in previous experiments [191 ] . 



VII. UPPER LIMITS FOR r/ c (25) PRODUCTION 

We search for the rj c (2S) in the invariant-mass distributions of the four- meson final states. 
Since the 7r + 7r~7r + 7r~ final-state sample has a prominent ip{2S) peak just above the r] c (2S) 
mass region, we veto events when the invariant mass of any 7t + 7t~ combination falls within 
0.1 GeV/c 2 of the nominal J/ip mass. (The two tracks from the J/ip decay are misidentified 
leptons.) 

We obtain upper limits on rj c (2S) yields from the four-meson invariant- mass distributions 
(Fig. 12). The mass and width of the r) c (2S) are not very precisely determined so far. We 
consider a wide range for the mass and width, 3.62 GeV/c 2 < M(r] r (2S)) < 3.67 GeV/c 2 
and 10 MeV< T(j] c (2S)) < 40 MeV, according to previous measurements [l5, 16|, 1^, 25, 26 



Fits similar to those made in Sect. III. A are applied for the rj c (2S) using a Breit-Wigner 
function with the mass resolution fixed to 9 MeV and a second-order polynomial for the 
background component, as well as a Gaussian function for the Xc2 peak. The upper limits 
for the yields are 340, 164 and 55 events at 90% C.L., for the Air, 2K2n and AK final-states, 
respectively. The curves corresponding to these upper limits are shown by the solid curves 
in Fig. 12. We calculate the upper limits using x 2 values from fits. We assume various 
non-negative values of the r) c (2S) yield, and, for each case, we obtain the x 2 from the best 
fit while floating all the other parameters. We define the 90% C.L. upper limit of the yield as 
one whose x 2 is larger by (1.64) 2 than the minimum x 2 derived in the different assumptions 
of the yield. An application of the Feldman- Cousins method (Table X in Ref. [20]) to the 
yields obtained from the fits gives the values very close to the yield upper limits, 338, 163 
and 55 events for the Ait, 2K2it and AK modes, respectively. 

We calculate the efficiency assuming a phase space distribution. The upper limit for 
Q = r 77 £> for each decay mode of the r] c (2S) is shown in Table IIHI We find that the ra- 
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TABLE III: The results in this paper for Q = T^B and comparisons with previous measurements 
(shown as "direct" and "indirect" below, see text) 0, 0, The results with inequalities cor- 
respond to 90%-C.L. upper limits, where we include the systematic errors by shifting the upper 
limits upwards by la(sys). The "isospin -BF" factors are used to obtain the Q results from the 
experimental measurements. 



Process 


Isospin •Si' 1 


This paper 


Direct 


Indirect 




factor 


G (eV) 


G (eV) 


G (eV) 


Tj c — > tt+tt— tt+tt— 


1.0 


40.7 ± 3.7 ± 5.3 


180 ± 70 ± 20 


83 ± 24 


rj c — > K+ K~ ir+ir— 


1.0 


25.7 ± 3.2 ± 4.9 


210 ± 70 


102 ± 30 


7j c -► K+K-K+K- 

lis 


1.0 


5.6 ± 1.1 ± 1.6 


280 ± 70 


11 ± 5 


Tin > 00 


0.333 


< 39 




130 ± 43 


fir —> f 9 f 9 


0.320 


69 ± 17 ± 12 




74 ±36 


r/ c -> K*K* 


0.333 


32.4 ±4.2 ±5.8 




66 ± 22 


Tin — > f") fo 
' IC J Z J 2 


0.251 


49 ± 9 ± 13 






Tin — > d><b 

IC ^r^r 


0.243 


6.8 ± 1.2 ± 1.3 




19 ± 5 


Yr-n — ► 7r + 7r~7r + 7r~ 

/V6U 


1.0 


44.7 ±3.6 ±4.9 


75 ± 13 ± 8 


69 ± 13 


XcO —> K+K-TT+TT- 


1.0 


38.8 ±3.7 ±4.7 




53 ± 12 


y K+K-K+K- 


1.0 


7.9 ± 1.3 ± 1.1 




7.8 ± 1.6 


XcO — ¥ K*°K~ir + or c.c. 


0.667 


16.7 ±6.1 ±3.0 




34 ± 13 


XcO ~^ PP 


0.333 


< 12 






Xco -> K*K* 


0.333 


< 18 




5.1 ± 1.9 


Xco -> 4>4> 


0.243 


2.3 ±0.9 ±0.4 


— 


2.7 ±0.8 


Xc2 TT+TT-TT+TT- 


1.0 


5.01 ±0.44 ±0.55 


6.4 ± 1.8 ±0.8 


7.2 ± 1.2 


Xc2 —> K+K~ir+TT~ 


1.0 


4.42 ±0.42 ±0.53 




5.8 ±2.1 


Xc2 -> K+K-K+K- 


1.0 


1.10 ±0.21 ±0.15 




1.03 ±0.18 


Xc2 -> p°7r + 7r _ 


1.0 


3.2 ±1.9 ±0.5 




3.9 ±2.3 


Xc2 ~> PP 


0.333 


< 7.8 








0.333 


2.4 ±0.5 ±0.8 




2.2 ±0.5 


Xc2 ~> # 


0.243 


0.58 ±0.18 ±0.16 




1.0 ±0.3 


r] c (2S) — ► 7r + 7r _ 7r + 7r~ 


1.0 


< 6.5 






?7 C (25) ^r+K-TT+vr- 


1.0 


< 5.0 






r/ c (25) -> K+K-K+K- 


1.0 


< 2.9 







tios Q(rj c (2S) —>■ four mesons )/Q{Vc{ IS) four mesons) are much smaller than unity in all 



three processes, as previously found in another decay mode, r} c (lS, 2S) — ► KgK+n^- jl5 
Although these results are still marginally consistent with the theoretical expectations, 
r 77 (?7 c (2S'))/r 77 (?7 c (lS')) ~ 0.3 and B(r) c (2S) — > some hadronic final state) / B(r] c (lS) — > 
the same final state) ~ 1, they give significant constraints on theoretical models [1-3], in- 
cluding calculations based on relativistic qq production. 
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FIG. 12: The invariant-mass distribution in the vicinities of rj c (2S) mass in the final states, (a)47r, 
(b)2K2n and (c)4if. The histogram and plots in (a) are the distributions before and after the 
rejection of ip(2S) — > J/-07r + 7r~, respectively. The dotted and solid curves show the best fit to 
the experimental data and the fit with the rj c (2S) yield set to the obtained 90%-C.L. upper limits, 
respectively. 

VIII. CONCLUSION 

The production of the T] C (1S), Xco and Xc2 charmonium states in two-photon colli- 
sions has been observed in all of the four-meson final states, 7r + 7r~7r + 7r~, K + K~7t + tt~ and 
K + K~K + K~ . We used data samples with one or two orders of magnitude larger statistics 
than previous measurements. No clear signature for the r] c (2S) is found in any of decay 
processes, and we obtain the upper limits for the products of its two-photon decay width 
and the branching fractions. We have studied resonant substructures in these four-meson 
final states. 

For the first time Xcj signals produced in two-photon collisions are observed in 
the K + K~7t + tt~ or K + K~K + K~ final states. We also find a new decay mode, 
r] c — > / 2 (1270)/2(1525). We have obtained products of the two-photon decay width and 
branching fractions for various decays of charmonium states. The present results for r] c 
are systematically smaller than the derived values from the world averages of previous 
measurements. 
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